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Abstract 

We present an analytical formula for the Bose-Einstein correlations (BEC) which 
includes effects of both Coulomb and strong final state interactions (FSI). It was 
obtained by using Coulomb wave function together with the scattering partial wave 
amplitude of the strong interactions describing data on the s-wave phase shift. We 
have proved numerically that this method is equivalent to solving Schrodinger equa- 
tion with Coulomb and the s-wave strong interaction potentials. As an application 
we have analysed, using our formula which includes the degree of coherence and the 
long range correlation, the data for e + e~ annihilations. 

We have found that the degree of coherence present in our formula approaches 
approximately unity whereas the long range correlation parameter becomes approx- 
imately zero. These results suggest that the physical meanings of the fractional 
degree of coherence and the long range correlation observed in various other analy- 
ses can most probably be attributed to FSI. 
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1 Introduction 



Very recently several analytical expressions for the final state interaction (FSI) corrections 
to Bose-Einstein correlations (BEC) due to Coulomb interaction have been presented in 
|l|, |2|] (see also ||, U). On the other hand it is also well known that there exists FSI due 
to the strong interactions between identical pions 0, ||, [jj. In order to treat both the 
FSI of the Coulomb type and those emerging from the strong interactions at the same 
time, one has to solve Schrodinger equation with the corresponding Coulomb and strong 
potentials. Such study has been made by Pratt et al. || x . However, in approach of 
H it is difficult to elucidate the physical meaning of the degree of coherence parameter 
introduced in various analyses of data on the BEC (see below). 



In this paper we consider this problem using different approach from that of ||. 
Our aim will be to obtain an analytical formula for the BEC which would include 
also FSI caused by both Coulomb and strong interactions. Comparing our results with 
those of || we prove numerically that both methods are equivalent. The data of the 
BEC in e+e" annihilation by TPC, AMY, OPAL ALEPH and DELPHI Collaborations 
p0| , |2"T| , p2| , p3| , p4| are then analysed using our formula. Finally concluding remarks are 
presented. 



2 Phase shift and wave function of identical tt-tt scat- 
tering 

Let us first consider the Coulomb wave function of the identical n-n scattering with 



momenta pi and P2 [12, 13 



*c(k,r) = L '(1 + ir))e-™ 1/2 ! e ik ' r 'F(-ir); 1; ikr (1 -cos 0)), (1) 



1 Study exploiting WKB method is presented in ||. See also (l(J and (IT), where pp correlations were 
studied. 
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(2krye lkr F(l + l + ir];2l + 2; -2ikr). 



(2/ + 1)! 

In numerical calculations of the hypergeometric function appearing in (1) wild (oscil- 
lating) behavior emerges when kr becomes large. To avoid the computational problem 
in this asymptotic region, we use there the following asymptotic formula of the Coulomb 
wave function: 



\l/c ym (k, r) = exp{i(kz + i]\n(k(r — z)))} 1 + 
, ffa\ ex v{ i ( kr ~r]\n(2kr))} 

+n u ) : > 



rj- 



ik{r — z) 



(2) 



where z = rcosd and the scattering amplitude f{9) is given as follows: 

T] 1 



f(0) 



exp{—2ir/ In sin(0/2) + 2i argT(l + if])}. 



2ksm 2 (9/2) 

The asymptotic formula for an s-wave component of the Coulomb wave function is ob- 



tained in the similar way |T3 



•o)( k : r ) = exp{i(kr -r]\n(2kr) + 2r?o)}^ 



+ exp{— i(kr — rj \n(2kr))}- 



-2ikr 



1 + 
1 + 



ir](l + irf) 
2ikr 

— irf) 
2%kr 



(3) 



As the next step let us consider the strong interaction provided in terms of the phase 
shift in the identical tt-tt scattering. The data for the s-wave (J = 2) phase shift reported 
in [14] -[19] can be described phenomenologically by the following formula ||: 



6 m 



1 



2 I 1 + 0.5Q 2 



(4) 



3 



where parameter ao ( -1.5 < ao < -0.7 (GeV 1 ) ) denotes the scattering length. 

Finally, using (4), the scattering s-wave function induced by strong interaction is 
expressed as: 

MKr) = f(0) eXP{i{kr - nH2kr))} , (5) 

r 

f°(0) = ^exp(2^ )(exp(2^f) - 1). 

It should be noticed here (cf. Fig. 1(a)) that in addition to (5) we need one more it-it 
scattering wave function, which would reasonably describe the strong interaction in the 
small kr region. Of course Coulomb potential affects both regions, cf. an "unknown box" 
in Fig. 1(a). To obtain a wave function suitable for that box we use the following assump- 
tion proposed by Bowler |J 0: We assume that the wave function given by (5) with a 
renormalization provided by the square root of the Gamow factor can be interpolated into 
the internal region. This is attributed to the normalization of the Coulomb wave function 
which is given by the following factor: ^jG(k) = (2irri / (exp(2irri) — l)) 1 ^ 2 . 

For the asymptotic region we use the following expression: 

* total (k, r) = vD^ sym (k, r ) + st (k, r) (6) 

whereas for the region described by the exact Coulomb wave function the following wave 
function is used instead: 

*totai(k, r) = * c (k, r) + ^G(AO0 st (k, r). (7) 

As seen in Fig. 1(b), there is smooth connection between both regions. The usefulness of 

this assumption will show up in Section 4. 

2 Cf. also @ where the plane wave and the data of the s-wave phase shift of identical ir-ir scattering 
were used. In the present calculation we have to use the Coulomb wave function instead of the plane 
wave. 
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3 Formulation of BEC 



To describe a pair of the identical bosons, we have to symmetrize the total wave function 
in the following way: 

1 



A 12 = ^[*c(k, r) + tfg(k, r) + $ st (k, r) + $f t (k, r)], 



(8) 



where superscript S denotes the symmetrization of the wave function. The function 
$ st (k, r) stands for the wave function induced by strong interactions. Assuming a source 
function p(r) we obtain the following expression for the BEC including the FSI: 



J c = E 



G(k) 

Ic + ^Cst + 4t, 

1 



m,n=0 



X 



m + n + 1 



n\m\ 



I R1 (2 + m + njA^rijAKm) 



n 
in 



(n + m)\ \ 
1 + Aic) + (E 2B + A EC 



m 
vq 



lest = 25R 



-(2A;)^exp (-i(r) + ^ 2) )) sin^ 2) £ ^(1 + n)A 2 (n, 0) 



n=0 



/s* = ^/Ri(0)sin 2 ^ 2) , 



(9) 

(10) 

(11) 
(12) 

(13) 



where 



E 2 b 
1 + A 1C 

A 1 (n) 
A 2 (n,l) 
Im(n) 



= J d 3 r p{r) e 



iQ r 



1 + An ■ 2r] / pr 2 dr ^ 



n A 2n+l 



f (2n+ l)!(2n + l)(2ra + 2)' 



r(n + ir?) (-2ik) n 
T(ii]) (n!) 2 ' 
r(n + j + 1 + i?7)r(2Z + 2) (-2ifc) n 
r(Z + 1 + ir])F(n + 21 + 2) n\ : 

47r y drr n p(r), 
4tt y drr n+ >(r). 



In this paper N( ±zL yN BG stands for the ratio of pairs of identical charged bosons in 
a single event to those from different events. Whenever the data are corrected by the 



Gamow factor the final formula should also be divided by the Gamow factor. In numerical 
computations we have to assume explicitly some forms of the source function. In present 
calculation we use the gaussian source function: p{r) = (^i^) 3 exp(^). Therefore the 
Fourier transform of the source function is given as follows: 

E 2B = exp(-/5 2 g 2 /2). 



4 Comparisons of our results with those of Pratt et 
al. 

Authors of H have presented their results for the BEC with interaction ranges (3 = 2 fm 
and (3 = 20 fm (corresponding to R in their notation) by solving the following Schrodinger 
equation: 

r 79 r / T , + \ "1 

li(x) = 0, (14) 



e 



dx z x A x 



where \i = m/2, x = //r, K 2 = 2E/fi, e = 2a and U(x) = (2Vo/m p x) exp(— m p x/p) for 
the s-wave. The values of parameters: Vq = 2.6 GeV and m p = 0.77 GeV were used. In 
Fig. 2 we compare our results with theirs using the same values of parameters {(3 = 2 fm 
and 13 = 20 fm). Because 7c with f3 = 20 fm shows wild oscillations near 50 MeV/c, which 
are due to the series expansion of the confluent hypergeometric function (n + m = 75), 
in our calculation we have used a method of seamless fitting introduced in ||. If we use 
a different constraint, n + m = 50, we observe a sharp decreasing near 42 MeV/c, all 
this depends on the parameter (3. The origin of this phenomenon can be attributed to 
the mathematical property of the confluent hypergeometric functions and the ability of 
computers. Therefore in the asymptotic region (2) should be used. 

As seen in Fig. 2 our results and solutions of the Schrodinger equation || are numeri- 
cally equivalent to each other except for the behavior near Q = a few MeV which is due to 
the logarithmic term in (5). It is therefore confirmed that from their potential (Vq = 2.6 
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GeV) we obtain in the Born approximation the scattering length ao 



= -0.8 ~- 0.6 GeV" 1 . 



In the actual analyses we have to introduce a cutoff parameter in the small Q region. 
Since in many cases there are no data, or available data have large error bars due to 
limits of momentum resolutions, analyses of data do not critically depend of this cutoff 
parameter. 



5 Analyses of data in e + e annihilation 

As stressed in Sec.l, we want to elucidate the physical meaning of the degree of coherence. 
The parameter A describing it should be therefore introduced into (9) in the usual way. 
Moreover, notice that two more parameters: the additional normalization factor c and 
the long range correlation parameter 7 are also introduced by hand. Our final formula 
containing all these parameters is thus given as: 



jV(±±)/jV BG (Q = 2k) = c(l + A 1C + A EC + I Cst + J st ) 

E2B 



x 



1 + A- 



I + 7Q). (15) 



1 + A 1C + A EC + / Cst + / St . 
It should be noted that the normalization c and an effective degree of coherence, i.e., the 
denominator of the ratio E 2 b/(^ + Aic + Aec + -^Cst + 1st), are related to each other. 



For the sake of reference we use in our analyses also the conventional formula (i.e., the 
standard formula without corrections due to the FSI): 

N (±±:Standard) /iV BG (Q = 2k) = c [1 + \E 2B ] (1 + 7 Q). (16) 



We apply our formulae to data for e + e~ annihilation provided by P0|- |P3 |. Results of our 

analyses performed by means of (15) and (16) are shown in Fig. 3 and Table I. f] As seen in 

Table I, estimated values of the degree of coherence parameter A are systematically larger 

(approaching unity) when (15) is used than those obtained by the standard formula (16). 

3 It is difficult to treat ao as a free parameter in the CERN MINUIT program, due to the limited 
ability of our computer. 



On the other hand, estimated values of the long range correlation parameter 7 approach 
approximately zero when (15) is used (except for the result found in present analysis of 
data by AMY collaboration). 



Finally please notice that reported values in |2(| -[24| are obtained by the following 
formula: 



^(d=h:Empiricai.)/jv BG (g = 2k) = c \l + X exp(-Q 2 R 2 )] (1 + jQ) 



(17) 



To compare our estimated values with various reported ones, we have to use a relation, 
(3 = \/2R, because of the different method of integrations. In Table I we show therefore 
the corrected values instead of the values reported by various collaborations. The /3's 
obtained in the analysis by (15) are systematically smaller than the corrected values 



taken from [j30| -|2|] and the estimated ones obtained in the analysis by the standard 
formula (16). 



6 Concluding remarks 

We obtain analytic formula (9) for BEC including the Coulombic and strong FSI. It is 
(numerically) confirmed that our method is equivalent to solving the Schrodinger equation 
(14). 

Combining the seamless fitting method and the CERN MINUIT program in (15) 
we have analysed data for BEC in e + e~ annihilations. Our results are significantly dif- 
ferent from those obtained by the standard formula (16). Namely, it is found that the 
degree of coherence parameter A and the long range correlation parameter 7 approaches 
approximately unity and zero, respectively (see Table I). Therefore, we conclude that the 
physical meanings of the fractional degree of coherence parameter A and the long range 
correlation parameter 7 as obtained by the standard formula (16) should be attributed 
to the FSI. Moreover, the values of the source size parameter reported by various collab- 
orations (after using the relation: (3 = y/2R) and obtained by the standard formula (16) 
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are systematically larger than values estimated by (15). 
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Figure Captions 



Fig. 1. (a) Interrelation between the two wave functions in external and internal 
regions, (b) Real part of total wave function with Q = 100 MeV/c and 9 = tt/2. See (7). 
Near 250 MeV/c there is a connecting point. 

Fig. 2. Comparisons of our results (solid line) and those of Pratt et al., || (dashed 
line). The scattering length ao = —0.6 GeV -1 is used. This value approximately corre- 
sponds to Vo = 2.6 GeV in the Born approximation, (a) /3 — 2 fm and (b) j3 = 20 fm. 
Behaviors in small Q region (see inserts) are attributed to the logarithmic term in (5). 



Fig. 3. (a) Analysis of data of TPC Collaboration |2(| by (15) and (16). 

(b) The same as (a) but for data of AMY Collaboration (The point at the smallest 
value of Q is neglected in analysis by means of (15).). 

(c) The same as (a) but data of OPAL Collaboration [^| . 

(d) The same as (a) but data of ALEPH Collaboration ||23|| . 

(e) The same as (a) but data of DELPHI Collaboration [El . 



Table Caption 

Table I. Analyses of data of the BEC by TPC, AMY, OPAL, ALEPH and DELPHI Col- 
laborations; ao=-1.00 GeV -1 . The source size parameters obtained by all collaborations 
are corrected by (3 = \[2R. There is no significant difference between those corrected 
values (denoted by (*)) and values obtained in present analysis. The AMY collaboration 
has used the fitting function R mix (Q) = C(l + f v {Q)X exp(-/? 2 g 2 /2))(l + ^Q),U{Q) = 
0.719 - 0.070Q + 0.056Q 2 - 0.020Q 3 . 



12 



P[fm] 



A 



7 



c x 2 / NDF 



TPC 

[20]: 
Our analyses: 



0.92 ±0.06* 0.61 ±0.05 

(15) 0.737 ±0.050 1.097 ±0.042 

(16) 0.912 ±0.062 0.611 ±0.054 



-0.000 ±0.020 1.002 ±0.022 44.2/35 
0.083 ±0.025 0.881 ±0.023 41.0/35 



AMY 
[21]: 
Our analyses: 



0.823 ±0.088* 0.392 ±0.041 

(15) 0.460 ±0.039 0.947 ±0.033 

(16) 0.854 ±0.088 0.286 ± 0.030 



0.033 ±0.041 0.935 ±0.016 90.2/93 
-0.078 ±0.012 1.203 ±0.030 106.8/95 
0.030 ±0.013 0.956 ±0.016 94.1/96 



OPAL 

[22]: 
Our analyses: 



1.124 ±0.021* 0.846 ±0.025 

(15) 1.090 ±0.035 1.043 ±0.025 

(16) 1.339 ±0.035 0.713 ± 0.036 



336/73 

0.003 ±0.004 0.991 ±0.005 124.4/74 
0.040 ±0.004 0.936 ±0.004 118.7/74 



ALEPH 



1.14 ±0.06* 0.48 ±0.03 
Our analyses: (15) 0.917 ±0.037 1.070 ±0.024 
(16) 1.128 ±0.037 0.630 ±0.030 



0.02 ±0.01 0.97 ±0.01 77/70 
-0.016 ±0.008 1.032 ±0.011 89.0/69 
0.024 ±0.009 0.964 ±0.011 87.3/69 



DELPHI 

[24]: 
Our analyses: 



1.16 ±0.04* 0.45 ±0.02 

(15) 0.871 ±0.027 0.946 ±0.016 

(16) 1.170 ±0.039 0.451 ±0.020 



89/73 

-0.001 ±0.006 1.020 ±0.008 90.3/73 
0.033 ±0.007 0.963 ±0.008 89.1/73 
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